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ABSTRACT 


This Addendum describes the continuation of niimerical studies 
dealing with the dynamics of minor constituents in the earth’s thermo- 
sphere. Previous investigations of the distribution of thermospheric 
helium were repeated, making use of the Jacchia (19’7V) model in place 
of CIRA (1972) for specifying the background gas temperature and pres- 
sure. Changes in the global transport of helium under solstice conditions 
caused by a small increase in the latitudes at which the background gas 
pressure extrema occur lead to much better agreement of the model pre- 
dictions with data taken by the mass spectrometers on board the ESRO-4 
and OGO-6 satellites. Development of a three -component (N^j 
O ) model of the thermosphere, initiated under funding by Aerospace 
Corporate Programs for Research and Experimentation, wa s continued. 

The model was applied to a study of the global distributions of these three 
atmospheric gases at both equinox and solstice, with emphasis on the winter 
enhancement of atomic oxygen in the lower thermosphere. Comparison of 
the results with measurements taken by the ESRO-4 mass spectrometer 
indicates that the distribution of atomic oxygen can be generally understood 
as being a result of global transport by winds. Oxygen chemistry, not in- 
cluded in the calculations to date, could modify these results, but, because 
of the time scales involved, would not change the basic conclusions. 
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I. 


INTRODUCTION 


During the first year of this contract, calculations dealing with 
the global distributions of the minor gases He and Ar were carried out. 
The calculations made use of a numerical model of the thermosphere in 
which the influence of a background gas on the dynamics of a minor gas 
is considered. The model covers the range 90-500 km, and the back- 
ground gas density, mean molecular weight and temperature were taken 
from the CIRA (1972) empirical model. Wind fields in the background 
and minor gases are calculated using the equations of conservation of 
momentum, including the effects of. anisotropic ion drag, viscosity and 
the Coriolis force, as well as pressure gradients derived from the den- 
sity and temperature gradients. 

Although the previously described model calculations reproduced 

the general observed features of the "winter helium biilge, " including the 

decrease of its amplitude with increasing solar activity and the latitudes 

and local times of the helium density extrema, quantitative comparison 

of the results of the model calculations with data taken by satellite -borne 

mass spectrometers indicated that the computed amplitude of the winter 

helium bulge was somewhat underestimated. For example, for 
-22 2 

Ffo y = 92xl0 W/m Hz and Kp = 0, our model predicted a helium 
bulge ratio (ratio of maximiim to minimum helium density at a given 
altitude) of ~ 28 at 275 km altitude, as opposed to a value of — 42 
measured by the ESRO-4 satellite (Keating et al. , 1976; von Zahn et al. , 
1977). Similarly, for Fj^q j = 140 x 10 W/m^ Hz, the theoretical 

model predicted a helium bulge ratio of 12 at 400 km, whereas the 

OGO-6 mass spectrometer gave a value 20 (Hedin et al. , 1974). 
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Since the primary mechanism responsible for the winter helium 
bulge is transport of helium from the summer to the winter hemisphere 
by global-scale winds in the background gas, the underestimate of the 
amplitude of the winter helium bulge is attributed to an underestimate 
of the speed of the meridional wind flowing from the summer to the 
winter hemisphere. This wind speed would be increased if the latitudes 
at which the extrema of the background gas pressure occur were larger 
than those given by the CIRA (1972) model. The CIRA (1972) model 
places the pressure maxim\xm at the sub solar latitude ( 23. 5° tmder 
summer solstice conditions ) and the pressure minimum at the anti- 
subsolar latitude. That these are the correct locations for the pres'sure 
extrema is far from clear. Data taken by the satellite -borne mass 
spectrometers mentioned earlier indicate that the solstitial temperature 
extrema occur at considerably higher latitudes; the locations of the density 
extrema depend on the molecular weight of the gas in question, and the 
locations of the total pressure extrema generally differ from those 
given by the CIRA model. 

Recently, Jacchia (1977) has presented a new empirical model of 
thermospheric temperature, density and composition. This model in- 
cludes changes in composition observed by the mass spectrometers on 
board the OGO-6 and ESRO-4 satellites, while keeping the total density 
profiles anchored to satellite drag data. The model has a formulation 
-.- ■basically similar to that of CIRA (1972), but its detailed model predic- 
tions are quite different in several respects. Of most relevance here is 
the fact that, even under = 0 conditions, the exospheric temperature 
extrema at solstice occur at ~ ±45° latitude, considerably poleward of 
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their positions in CIRA (1972). Thus, even though the total density 
extrema occur at latitudes near the subsolar and anti-subsolar latitudes, 
the total pressure extrema occur somewhat poleward of the correspond- 
ing positions in CIRA (1972). Thus, one should expect that use of the 
Jacchia (1977) model, in place of the CIRA model, to specify the 
temperature and background gas pressure would lead to an increase 
in the helium bulge ratio, all other factors remaining unchanged. 

This feature of the Jacchia (1977) model has a further relevant 

effect. Exospheric flow of helium tends to reduce the amplitude of the 

winter helium bulge. The amount of exospheric flow depends on the 

2 5/2 

value of the quantity Q = V (nT ) evaluated at the base of the exo- 

X 

sphere; here n is the helium number density, T is the temperature, 

2 

and V is the horizontal Laplacian operator. Since the maximum of n 

X 

occurs at high latitudes in the hemisphere in which T takes on its minimum 
value, Q is reduced by increasing the latitudes at which the extrema of 
T occur. Thus, use of the Jacchia (1977) model is expected to lead to an 
increase in the predicted amplitude of the winter helium bulge, both 
because it increases the flow of helium from the sxxmmer to the winter 
hemisphere in the thermosphere, and because it decreases the return 
flow of helium in the exosphere. 

Section II of this Addendum describes an investigation of the extent 
to which use of the Jacchia (1977) model can improve the agreement be- 
tween the predictions of our theoretical model and measurements of the 
amplitude cfthe winter helium bulge measured by the ESRO-4 and OGO-6 
mass spectrometers. 
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In addition to the inert gases, which are minor constituents of 
the thermospheire at all altitudes below ~ 500 km,, atomic oxygen is 
a minor species in the lower thermosphere (below ~ 140 km). That the 
distribution of atomic oxygen in the thermosphere is sensitive to pheno- 
mena other than photochemical processes and hydrostatic equilibrium 
was first suggested by King (1964) and Johnson (1964). They discussed 
the importance of horizontal transport of atomic oxygen from the summer 
to the winter hemisphere and how it can explain the observed decrease 
in the fraction of molecular ions in the F -region in winter. Several 
sources of data have confirmed this conclusion. Measurements using 
the satellites Explorers 19 and 39 near; opposite poles at 1000 km 
altitude were used by Keating et al. (1971) to infer that the maximum 
O concentration near 120 km altitude occurs at low latitudes in 
the winter hemisphere. They pointed out that this winter enhancement 
of O can result from the same type of meridional flow into the winter 
hemisphere that produces the winter helium bulge (Keating and Prior, 

1967, 1968). The fact that, at higher thermospheric altitudes, the 
atomic oxygen maximum occurs in the summer hemisphere was attri- 
buted to the fact that temperatures in the summer hemisphere exceed 
those in the winter. Barlier et al. (1971) and Alcayde et al. (1974) 
used incoherent scatter radar measurements and satellite drag data 

i 

to infer=,a winter maximum of atomic oxygen in the lower thermosphere. 
Donahue et al. (1973) gave an exhaustive discussion of OGO-6 observations 
of the ‘5577 nightglow and their implications for the atomic oxygen dis- 
tribution in the 80-120 km altitude range. They demonstrated that, 
assuming that latitudinal variations in the eddy diffusion coefficient 
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are not large enough to be dominant, the large latitudinal variations in 
the density of O requires meridional flows of 10-50 m/sec at the 100 km 
level. Offermann (1974) presented a collection of density and composition 
data taken by instruments on various rockets; he showed that in the 120- 
200 km altitude range, the minim.um O density occurs ~ 30 days after 
summer solstice. He also demonstrated that at 150 km, the seasonal 
variation of the n(0)/n(0^) ratio is 2 and is thus larger than vari- 
ations correlated with solar flux, geomagnetic activity or local time. 

More recently, mass spectrometer data have been used to describe the 
global variation of atomic oxygen concentrations. (See for example, 

Hedin et al, , 1973; Newton et al. , 1975; Keating et al. , 1976; ktauersberger 
et al. , 1976;von Zahn et al. , 1977). These data have been incorporated 
into the Mass-Spectrometer Incoherent Scatter Model by Hedin et al. 

(1977). 

Theoretical studies dealing with the global distribution of atomic 
oxygen include the seasonal model of Mayr andVolland (1972), in which 
both helium and atomic oxygen winter bulges were found to be produced 
in the lower thermosphere by circulation. At higher altitudes, thermal 
expansion was shown to be more important than circulation for O, leading 
to a maximiim of O in the summer hemisphere above ~ 440 km altitude. 
Johnson and Gottleib (1973) and Johnson (1973) demonstrated that a 
meridional wind system that is required from thermal budget consider- 
ations is adequate to transport the needed quantity of O from the summer 
to the winter hemisphere. Johnson concluded that, although the same 
transport processes act on O as on He, the effect on O should be 


5 



smaller than that on He because the scale height of O is less than a 
factor of Z greater than the average atmospheric scale height in the 
lower thermosphere. Finally, Mayr and Harris (1977) presented a 
theoretical model of the diurnal variations in thermospheric tempera- 
ture, composition and winds. Using perturbation theory, they showed 
that thermospheric winds have a substantial effect on atomic oxygen 
below 200 km. They tend to shift the diurnal maximum of n(0) in the 
lower thermosphere to times earlier than would occur under the action 
of hydrostatic equilibrium alone. A gradual phase transition occurs, 
with maximum n(0) oc curing at ~ 6 LT at 150 km and at ~ 14 LT at 
300 km. Chemistry was shown to be of secondary importance in de- 
termining the diurnal variations of n(0). 

As part of our investigations this year, a study of the dynamics 
of atomic oxygen has been carried out. Because atomic oxygen becomes 
the major thermospheric constituent above ~ 150 km, the numerical 
model used in the study of He and Ar could not be used to study the 
distribution of atomic oxygen. An interactive three -component model, 
whose development was begun under funding by Aerospace Corporate 
Programs for Research and Investigation, was used to investigate the 
dynamics .of N 2 > and O . The mathematical formulation and imple- 

mentation of this model are described in Sections III and IV. 
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II. 


USE OF THE JACCHIA (1977) MODEL IN DETERMINING THE 
GLOBAL DISTRIBUTION OF THERMOSPHERIC HELIUM 


The complete formulation of the three-dimensional theoretical 

model was described in the original report (Straus et al. , 1977). The 

only modification made for the present calculations involved using 

the Jacchia (1977) model in place of the CIRA (1972) model for specifying 

the temperature and backgroimd gas pressure in the 90-500 km altitude 

range. Calculations were carried out for June solstice conditions with 

-22 2 

values of F^^ ^ = 92 and 140 x 10 W/m Hz; these values allow 

direct comparison to be made with the data taken by the ESRO-4 and 

OGO-6 mass spectrometers. In order to isolate the effects of changes 

in the background gas distribution, all calculations were carried out 

6 Z 

for K =0 and eddy diffusion coefficient K =6x10 cm / sec. In 
P 

describing the results of these calculations, we shall refer to the use 
of CIRA (1972) as "Model 1" and the use of Jacchia (1977) as "Model 2". 

Figure 1 shows the zonally averaged meridional velocity v^ 

across the equator for Models 1 and 2 for the two values of F^q 

(PDsitive Vq corresponds to southward flow. ) In ail cases, throughout 

most of the thermosphere the zonally averaged motion is directed 

toward the winter hemisphere; a small region of reversed flow occurs 

at low thermospheric altitude. Generally, v^ increases with Fj^q 

because of an increase in the ion density with F^q , Vg does not 

increase with Fjg ^ as rapidly as do- the pressure gradients. For 

, both values of F,« ^ , the values of v^ calculated using Model 2 exceed 

it)*/ u 

those from Model 1; however, the increase is no more than 10 m/sec, 
and values of this magnitude occ,ur only in the upper thermosphere, where 
wind-induced departures from hydrostatic equilibrium do not occur. 
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3 

Figure 2 shows contours of constant values of log^^ { > g/cm ) 

-22 2 

at 275 km at June solstice for Fjq ^ = 92 x 10 W/m Hz calculated 
using both Models 1 and 2. Both figures show isopycnics which are 
strikingly similar in their general features and shapes to those generated 
by von Zahn et al. (1977) from data taken by the ESRO-4 mass spectrom- 
eter. Comparison of these two figures and a more detailed tabulation of 
the results indicate the following; , 

1) Model 2 produces a larger ratio of maximum to minimum He 
density than does Model 1, although the global average He 
density is essentially unchanged. 

2) Model 2 produces helium density extrema at higher latitudes 
than does Model 1. 

3) Model 2 produces helium density extrema at very slightly 
later local times than does Model 1. 

All of these resTilts are to be expected. The helium bxilge ratio 
and the latitudes of the helium density extrema increase because of 
changes in the meridional flow. The local times of the helium density 
extrema increase because the Jacchia (1977) model places the tempera- 
ture maximum at -^17 hours local time (and total density maximum at 
1415 LT ), whereas the temperature maximum in the CIRA (1972) model 
occurs at 1415 LT, along with the. total density maximum. Thus, the 
total pressure exti*ema occur at slightly later local times in the Jacchia 
(1977) model than in the CIRA (1972) model. 

Figure 3 compares the calculated vertical distribution of the 

helium bulge ratio with those given by the ESRO-4 empirical model 

-?? ? 

for Fj^q .j, = 92 X 10 W/m Hz and the OGO-6 empirical model for 
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Fjq ^ = 140 X 10 W/in . Hz. In each of the curves representing 
theoretical results, wind-induced departures from, hydrostatic 
equilibritim cause the helium bulge ratio to rise rapidly from its 
value 1 at 90 km to a -maximum near 250 km. Above that altitude, 
hydrostatic equilibrium prevails because molecular diffusion effects 
dominate large-scale transport processes, and the bulge ratio 
decreases at a rate dependent on the temperature. Thus, the 
behavior below ~ 250 km is that which the empirical models attempt 
to represent by spatially -varying lower boundary conditions. For a 
fixed value of F^^ the bulge ratio calculated using Model 2 exceeds 
that resulting from Model 1 at all altitudes and leads to substantially 
better agreement with the empirical models. Comparison of the OGO-6 
model for F^^ .^ = 140 x 10 W/m Hz with the theoretical calculations 
indicates that Model 2 yields a helium bulge ratio within 'v- 15% of the 
OGO-6 value in the 400-450 km altitude range, where the OGO-6 data 
were taken. Similarly, the predictions of Model 2 are within 15% 
of the ESRO-4 model for F^^ ^ = 92 x lO"^^ W/m^ Hz at 275 km. 

In summary , we have shown that use of the Jacchia (1977) model 
to specify the temperature and- background gas pressure distributions 
in a three-dimensional dynamical model of thermospheric helium leads 
to better agreement with observational data than-does use of CIRA (1972) 
(Straus et al. , 1977). The calculated increases in the amplitude of the 
winter helium bulge and the latitudes at which the helium density extrema 
occur are primarily due to increased meridional transport of helium 
from the summer to the winter hemisphere by large-scale winds. How- 
ever, the fact that the winter helium bulge amplitude decreases with 
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increasing solar activity is. not due to a decrease in wind speed with 
solar activity, but rather to the increased effectiveness of exospheric 
transport, which carries helium back to the summer hemisphere. It 
is also worthy of note that the amplitude of the winter helium bulge 
appears to be an extremely sensitive fvinction of the background gas 
pressure distribution, so the observed substantial variations in the 
helium distribution with geomagnetic activity are quite understandable. 
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III. 


THE THREE -COMPONENT MODEL 


In order to study the dynamics of the thermosphere, the model 
treats the three main gases N^, O^, and O , coupled together by 
collisional and chemical processes. The conservation of mass and 
momentum equations are solved in a spherical shell with lower and 
upper boundaries at 90 and 500 km, respectively. As in our earlier 
work, the temperature field has been taken from CIRA (1972). Under 
steady- state conditions, -the equations of mass^ and momentum con- 
servation of the i^ species in a spherical ( r, 0, 0 ) coordinate 
system are: 


Q3p./S0 + V. (p.v.) = S. -L. 
1 1 1 


(1 


Q(5v./S0) + v..Vv. + 2Qx V. = 



+ K 





+ K 


- ( 1/p. ) Vp. + ( fi./p. ) V V. - V. j { V. - Vj ) 


-E 


V.. ( V - V ) 
IJ r-J 


(2 


Here, is the angular rotation rate of the earth, p is the density, 

V is the velocity vector, g is the acceleration of gravity, p is the 

pressure ( p^ = p^kT/m^ , where k is Boltzmann's constant, T is 

the temperature and m^ is the gas molecular weight ), p is the viscosity. 
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V- T is the neutral-ion collision frequency ( v. t =2.6x 10 M n. , - 

1 1 It 1 

where M is the gas molecular weight in AMU and n^ is the ion number 
density, taken from the empirical model of Ching and Chiu (19’73) )v ^ 
is the ion velocity, taken to be equal to { v.b ) b , where v is the 
mass -average neutral velocity and b is the unit vector in the direction 
of the geomagnetic field. 

The last term in equation 2 represents momentum transfer due 
to collisions between neutral species. We use 

V-. = kTn. /m. D.. N, (3 

ij J 1 U 

where N = E n, and D.. is the binary diffusion coefficient (Banks and 
k ■ U 

Kockarts, 1973). The gravity term allows a transition from complete 
mixing in the lower thermosphere ( K > > D. ) to diffusive equilibrium 
in the upper thermosphere ( D. > > K ). Here, 


D. 

1 
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All of the calculations to be described here were carried out with a 

6 2 

single value of K = 6 x 10 cm /sec. This circumstance is primarily 
due to the expense involved in the computations. Furthermore, as 
shown in our earlier work, a change in the globally uniform value of 
K does not have an important effect on the global variability of the 
thermospheric He density. Since the gases treated here have mole- 
cular weights much closer to the average than does helium, the re- 
sults should be even less sensitive to the value of K used, as long 

6 T 2 

as it lies in the range 10 - 10 cm /sec. 

The terms S. and L. in { 1 ) represent source and loss terms 
due to chemical proces ses (photodissociation, photoionization, recom- 
bination, etc. ) The only processes of possible importance in the pres- 
ent model involve oxygen chemistry. Atomic oxygen is produced by 
photodissociation of in the lowest part of the thermosphere by 

solar radiation with wavelength k < 2424A. Because the recombination 
of O to form is a three -body reaction, photochemical equilibrium 

is accomplished by eddy transport of atomic oxygen from its point of 
formation downward into the mesosphere, where recombination occurs. 
This is an extremely complex process, and only investigations using 
one -dimensional photochemical models have been carried out. Under- 
standing of this process has not progressed to a point where oxygen 
•chemistry can be included in the present study. It is fortunate that 
the inclusion of oxygen chemistry is not crucial in understanding the 
global distribution of O in the thermosphere. The time required for 
photochemical eq\iilibrium to occur at 100 km altitude is on the order of 
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1 month and increases with altitude (Banks and Kockarts, 1973). Thus 
dynamical processes with time scales less than this will dominate the 
photochemical processes. This is certainly the case for diurnal 
variations. As will be seen, the winter enhancement of atomic oxygen 
in the lower thermosphere can be \mderstood as being a manifestation of 
meridional winds, whose time scale for transport from one hemisphere 
to the other is on the order of i-4 days. 

Equations 1-3 are solved simultaneously for the 3 gases 
using the spectral method described by Creekmore et al. (1975). Be- 
cause of the large size of the system of equations involved, special 
procedures for solving the coupled nonlinear algebraic equations were 
devised, tested and used extensively in this investigation. 

The boundary conditions employed were that, at 90 km, all 
velocities vanish and the densities of the 3 gases are those given in 
CIRA (1972). At 500 km, the horizontal velocities v. are assumed 
to have no vertical variation ( i. e. , Sv./Sz = 0 at 500 km ). 
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rv, RESULTS AND DISCUSSION 

Computations as .described in the previous section were carried 

-22 • 2 

out for equinox conditions with = 92 x 10 W/m Hz and for 

solstice conditions with F^^ ^ = 92 x 10~^^W/m^ Hz and 140 x 10~^^ 

2 

W/m Hz. The smaller value of ^ was chosen to allow direct 

comparison with the ESRO-4 measurements, whereas the larger value 
corresponds to those of OGO-6. The horizontal and vertical winds thus 
calculated are similar to those discussed in our first report (Straus et 
al. , 1977). In this report, we shall describe the density distributions 
of the three major gases in some detail, with emphasis on wind-induced 
departures from hydrostatic equilibrium. 

Figures 4-9 show contours of the natural logarithm of the density 

3 

(gm/cm ) of the atmospheric constituents and O at several 

altitudes under spring equinox conditions for Fjq ^ = 92x 10"^^ W/m^ Hz. 
Figures 4-6 show results for hydrostatic conditions^ whereas Figures 
7-9 illustrate complete calculations, including the effect of wind -induced 
diffusion. Comparison of corresponding figures for the hydrostatic 
and non -hydro static calculation shows important effects in the lower 
thermosphere. Molecular oxygen and nitrogen show some departure 
from hydrostatic equilibrium below ^ 200 km altitude, the departure 
being larger for O^, than for N 2 * Since has a molecular weight 

larger than the average (which is close to that of N 2 ijti the lower 
thermosphere) this is to be expected. The effects of winds are large 
enough to move the position of. maximum n( ) and n( ) away 
from the equator in the lower thermosphere and to generally change the 
shapes of the contours, although 1:he mean density at a given altitude is 
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essentially unchanged. Furthermore, the local times of maximum 
density (phase) of and are later in the non-hydro static cases 

than under hydrostatic conditions. Above 200 km, the effects of 
vertical motions on the and densities become negligible, 

and hydrostatic equilibrium prevails. 

Figures 6 and 9 show that departures from hydrostatic equili- 
brium are more substantial in the case of atomic oxygen than for N, 
or O^. Below 200 km, the atomic oxygen density has a pronoxmced 
maximum at 13 hours local time on the equator and strong minima 
at high latitudes in both hemispheres. The region from 200-300 km 
displays a transition to hydrostatic equilibrium, -with the effects of non- 
hydro staticity at lower altitudes being still visible at 470 km, especially 
in the asymmetry of the diurnal variation. 

Data relevant to global variation of the atomic oxygen density at 
equinox have been collected by several satellite -borne mass spectro- 
meters. Hedin et al. (1977) have used data from several satellites to 
generate the Mass Spectrometer -Incoherent Scatter (MSIS) Model. ' 
Comparison of this model which calculated resvilts presented here for 
atomic oxygen shows reasonably good agreement, especially in light of 
the uncertainties in several relevant input parameters, lower spatial 
resolution in the theoretical model, and extrapolations in the MSIS 
model to altitudes below those at which data were taken. Due to the 
neglect of oxygen chemistry, the global mean n( O ) is generally lower 
in the theoretical model than in MSIS. However, the global variations 
of n( O ) at a fixed altitude are similar in both magnitude and geographical 
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distribution. For example, at 120 km altitude, both the theoretical 
model and MSIS show minimum n(0) at high latitudes in both hemis- 
pheres and a broad region of larger n{0) near, and centered on, the 
equator. At 140 km, the MSIS model shows a 50% variation in n(0), 
whereas the theoretical model shows a 40% variation. Departures 
from hydrostatic equilibrium are clear in MSIS at all altitudes below 
~ 400 km, with a broad transition region in the 200-400 km altitude 
range. 

Figures 10-15 compare the hydrostatic and non-hydro static 
distributions of N^, and O at several altitudes for June solstice 
conditions with = 92 x 10 W/m^ Hz. As in the equinox cal- 
culations, the and distributions show some effects of wind- 

induced departures from hydrostatic equilibrium. In the lower thermo- 
sphere, winds drive the positions of the extrema of n( N_ ) and n( 0_ ) 
to higher latitudes and somewhat later local times than those calculated 
under the hydrostatic assumption. This behavior is similar to that 
described by Straus (1977) for thermospheric argon. Furthermore, 
the amplitudes of the global variations of n( N, ) and n( 0, ) are in- 
creased in the full calculation, relative to the hydrostatic calculation. 

For example, at 120 km altitude, the ratio of maximum n( N 2 ) to 
minimum n( ) , to which we shall refer as R( » is increased 

from 1.17 to ~ 1.30 ; similarly, R{ ) at 120 km is- increased 
from ~1.18 to 1.4Z. The global mean values of n( ) and n( O ) 
are essentially unchanged. At 185 km, the corresponding values for 
R( ), are I .90 and 2,30, and those for R{ ) are 2.04 and 
2. 79. 
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Above 200 km, molectilar diffusion effects become larger, 
and the and distributions become more nearly hydrostatic. 

The contours shown in panels d, e and f of Figures 10, 11, 13 
and 14 show that hydrostatic equilibrium clearly holds above this 
altitude. Furthermore, the scale heights of and are small 

enough that these species are quite sensitive to global temperature 
variations in the upper thermosphere. Thus, a transition region exists 
in the 200-300 km altitude range, above which the global distributions 
of n( ) and n( ) strongly resemble that of the temperature. 

Atomic oxygen is greatly affected by transport by global-scale 
winds,leading to a winter bulge of O in the lower thermosphere. Com 
parison of Figures 12 and 15 indicates that winds drive atomic ' 
oxygen from the summer hemisphere, where the n( O ) maximum 
occurs under hydrostatic conditions, to the winter hemisphere. The 
calculations indicate that this winter enhancement of O exists even 
at the lowest level in the model, 91 km, where the maximtim n( p ) 
occurs at - 30° latitude, 8 hours local time.. The global maximum 
n( O ) at June solstice at 120 km altitude occurs at ^10 hours local 
time, - 30 latitude. With increasing altitude, the maximum value 
of n( O ) at a given altitude moves towards the summer hemisphere 
and later local time. It reaches the equator at ~300 km altitude, 
where the local time of maximum n{ O ) is o- 13 hours. The slow 
transition with altitude of the location of maximum n( O ) is consistent 
with its relatively long scale height; atomic oxygen is less sensitive to 
global temperature variations than are and , 0 ^. (In the case 
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of the very light gas He, the transition effectively does not occur at 
all, and the maximum n { He ) occurs at' very high winter latitudes 
in the early morning hours. ) The amplitude of the winter enhancement 
of O , measured by R( O ), is shown in Table 1 as a function of alti- 
tude, for calculations with ^ = 92 and 140 x lO"^^ W/m^ Hz. 

(The contour plots for the latter value of y are given in Figures 

16-18. ) This table indicates that R( O ) below ^200 km increases 
slightly with F^^q y , but the increase is not significant. At higher 
altitudes, the ratio is less meaningful as a measure of the winter 
atomic oxygen bulge, since the extrema of n( O ) occur near the 
equator, and the winter bulge of O does not really exist. 

The results of these calculations can be compared directly 
with data taken by satellite -borne mass spectrometers. Von Zahn 
et al. (1977) give contour plots of the global distribution of N, and 
O at 275 km altitude iinder June solstice conditions with y = 

92x10 W/m Hz. Comparison of these figures with figures 13 d 
and 15 d shows very good agreement. For N^, the local times of 
the extrema of n( ) are faithfully predicted by the theoretical 
calculations. The latitudes of the extrema are somewhat closer to the 
equator in the theoretical model than in the empirical model; because 
is sensitive to the temperature distribution, this disagreement is 
probably due to the position of the temperature extrema in the CIRA ' 
(1972) model. • In the case of O , the agreement is quite excellent, 
considering that the small-scale features shown in the figure of von 
Zahn et al. cannot be resolved by the theoretical model. The local 
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TABLE 1 

Ratio of Global Maximum to Minimum n( O ) , R{ O ) 


Altitude ( km ) 

105 

119 

143 

185 

260 


R( O ), Fjq = 92 

1 . 22 
1.41 
1. 64 
1. 83 
1.91 


1 . 22 
1.46 
1. 73 
1. 89 • 
1, 73 


= 140 
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times of the extrema of n( O ) are well represented, and the lati- 
tudes of the extrema of n( O ) agree well. The value of R{ O ) is 
~2. 0 in both the theoretical and empirical models. 

Similar agreement exists between the theoretical model and 

the OGO-6 empirical model {Hedin et al. , 1974). Hedin et al. show 

contours of n( O ) at 460 km under June solstice conditions with 
-22 2 

Fjq ^ = 140 X 10 W/m Hz. Comparison of the plot with Figure 
18 f ( 470 km ) shows that the theoretical model predicts the correct 
location of the n( O ) extrema near the equator. Furthermore, the 
values of R( O } are very nearly equal, 2. 7 for the OGO-6 model 
and 2. 9 for the theoretical model. 

As mentioned earlier, the MSIS empirical model (Hedin et al. , 
1977) combines density data taken by several satellite -borne mass 
spectrometers. The general features of the distribution of atomic 
oxygen in the MSIS model are similar to those-of the ESRO-4 and 
OGO-6 models. The MSIS model, like our theoretical model, shows 
that the winter enhancement of atomic oxygen extends to very low 
thermospheric altitudes. The maximum n( O ) occurs at middle 
latitudes in the winter hemisphere, while the minimiim n( O ) occurs 
at high summer latitudes. In the 200-400 km region, the position of 
maximum n( O ) shifts towards the summer hemisphere, in response 
to the high temperature there. Finally, as in the theoretical model, 

R( O ) is relatively insensitive to ^ , increasing only slightly as 

Fj^q is increased from 90 to 150 x lO"^^ W/m^ Hz. (The global 
average value of n( O ), of course, increases with F^^ ^ .) 
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Comparisons of the results of the theoretical model with data 
that have not been included in global empirical models may also be 
made. There is a large amount of this type of data from both rockets 
and satellites, but most of it deals with absolute measurement of 
composition or ratios of the density of various species at a particular 
geographical locale, time of day, day of year, etc. The sources of 
data directly relevant to the present investigation, i. e. , the seasonal 
variation of the atomic oxygen distribution, are, in fact, quite limited. 
Offermann (1974), in a discussion of composition data obtained from 
rocket experiments in the altitude range 120-200 km, finds that there 
is a strong winter enhancement, with R( O ) ~ 1. 9 at 200 km. This is 
in good agreement ■with, the result of the theoretical model, as given 
in Table 1. Finally, Mauersberger et al. (1976) have discussed data 
taken at high northern latitudes in the winter and summer of 1974 by 
the open source neutral mass spectrometer on the AE-C satellite. They 
indicate that the ratio of n( O ) at ^ 60° latitude in early February to 
that in late June is 1. 5 at 300 km, eventually crossing 1.0 at 
^ 375 km. The results of the theoretical model are in resonably good 
agreement with these findings. 
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V. 


CONCLUSIONS 


In this report, we have described the continuation of theoretical 
studies of minor constituents in the thermosphere. Use of the Jacchia 
(1977) model for specifying the background gas pressure in our two- 
component model dealing with the distribution of thermospheric helium 
leads to considerably better agreement of the model predictions with 
observational data than we had obtained in previous calculations using 
CIRA (1972)., A new model of N^^ and O was applied to a study of 
the global distributions of these three gases, with particular emphasis 
on the seasonal variation of atomic oxygen. The model correctly pre- 
dicts a winter enhancement of atomic oxygen in the lower thermosphere, 
with the position of maximixm n(0) gradually shifting towards the 
summer hemisphere in the, thermopause region. Comparison of the model 
results with observational data shows generally good quantitative agree- 
ment, especially in light of the relatively low spatial resolution of the 
model. 
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FIGURE CAPTIONS 


Figure 1. 


Figure Z. 


Figure 3. 


Figure 4. 


Zonally averaged meridional wind speed at the equator 
calculated for F^q ,^ = 92 and 140 x 10"^^ W/m^ Hz 
using Models 1 and 2. 

3 

Contours of constant values of log^^Q ( , g/cm ) at 

275 km altitude at June solstice for F^q ^ = 92 x 10 
W/m^ Hz. 

a. Model 1 

b. Model 2 

Comparison of the He bulge ratio as a function of alti- 
tude as predicted by Models 1 and 2 for two values of 
Ffo y.with the empirical models based on data from the 
ESRO-4 and OGO-6 mass spectrometers. 

3 

Contours of Sm. ( , g/cm ) at spring equinox for 

hydrostatic conditions with F^^q ^ = 92 x 10"^^W/m^Hz. 


a. 

119 km 

b. 

143 km 

c. 

185 km 

d. 

260 km 

e. 

365 1cm 

f. 

470 km 
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3 

Figure 5. Contours of ^ , g/cm ) a.^ spring equinox for 

-22 ? 

hydrostatic conditions with ^=92x10 W'/m Hz. 

a. 119 km 

b. 143 km 

c. 185 km 

d. 260 km 

e. 36'5 km 

f. 470 km 

3 

Figure 6. Contours of Sn ( , g/cm ) at spring equinox for 

hydrostatic conditions with F^^q ^ = 92 x 10~^^/m^Hz. 

a. , 119 km 

b . 143 km 

c. 185 km 

d. 260 km 

e . 365 km 

f. 470 km 

3 

Figure 7. Contours of Sn ( p^^ g/cm ) at spring equinox for 

In ^ j 

^ -22 2 
non-hydro static conditions with ^ = 92 x 10 W/m Hz. 

a. 119 km 

b. 143 km‘ 

c . 185 km 

d. 260 km 

e. 365 km 

f. 470 km 
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3 

Contours of Sm. ( , g/cm ) at spring equinox for 

^ -Z2 2 

non-hydro static conditions with ^ = 92 x 10 W/m Hz. 

a. 119 km 

b. 143 km 

. c. 185 km 

d. 2 60 km 

e. 365 km 

f. 470 km 

3 

Contours of foz ( , g/cm ) at June solstice for hydrostatic 

-22 2 

conditions with .j = 92 x 10 Hz. 

a. 119 km 

b. 143 km 

c. 185 km 

d. 260 km 

e. 365 km 

f. 470 km 


3 

Figure 10. Contours of { p,, , g/cm ) at June solstice for hydrostatic 

^ -22 2 

conditions with ^=92x10 W/m Hz. 

a. 119 km 

b. 143 km 

c. 185 km 

d. 260 km 

e. 365 km 

f. 470 km 
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Figure 11. 


Contours of ( pQ , g/cm ) at June solstice for hydrostatic 
conditions with ^ = 92 x 10"^^W/m^ Hz. 

a. Il9 km 

b. 143 km 

c. 185 km 

d. 260 km 

e. 365 km 

f. 470 km 

3 

Figure 12. Contours of , g/cm ) at Jime solstice for hydrostatic 

conditions with F^^ ^ = 92 x 10“^^ W/m^ Hz. 

a. 119 km 

b. 143 km 

c. 185 km 

d. 260 km 

e. 365 km 

f. 470 km 

3 

Figure 13. Contours of S/n ( p^^ , g/cm ) at June solstice for non- 

^ -22 2 
hydrostatic conditions with ^ = 92 x 10 W/m Hz, 

a. 119 km 

b. 143 km 

c. 185 km 

d. 260 km 

e. 365 km 

f. 470 km 
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Figure 14. 


Figure 15. 


Figure 16. 


Contours of 0n { , g/cm ) at June solstice for 

^ -22 2 
non-hydro static conditions with Fjq y = 92 x 10 W/m Hz. 

a. 119 hm. 

b. 143 km. 

c. 185 km 

d. 260 km 

e. 365 km 

f. 470 km 

3 

Contours of £71 ( , g/cm ) at June solstice for 

-22 2 

non-hydro static conditions with F^^ y = 92xl0 W/m Hz. 


a. 119 km 

b. 143 km 

c. 185 km 

d. 260 km 

e. 365 km 

f. 470 km 


Contours of j0n ( , g/cm ) at June solstice for 

2 _22 Z 

non-hydro static conditions with F^q y = 140 x 10 W/m Hz. 


a. 119 km 
( 

b. 143 km 


c. 185 km 


d. 260 km 


e. 365 km 

f. 470 km 
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Figure 17. 


3 

Contours of ^ ( p , g/cm ) at June solstice for 
non-hydrostatic conditions with ^ = 140 x 10"^^W/m^Hz. 

a. 119 km 

b. 143 km 

c. 185 km 

d. 260 km 

e. 365 km 

f. 470 km 

Figure 18. Contours of , g/cm ) at June solstice for 

non-hydro static conditions with F^^^^ ^ - 140 x 10“^^/m^Hz. 

a. 119 km 

b. 143 km 

c . 185 km 

d. 260 km 

e. 365 km 

f. 470 km 
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-3.21317E.C 













latitude 





latitude 


OSPHERIC DE 


10.7 cn FLL 

I 


LOCAL TIME 


mBOL 

LN DCJ 

0. 

-2.52221E 

5. 

•2.52044E 

10. 

-2.51867E 

15. 

-2.51690C< 

20. 

-2.51514E< 

25. 

-2.51 337E 

30. 

-2.5H60C* 

35. 

-2.50983C 

40. 

-2.50806E 

45. 

-2.50629e 

50. 

-2.50452C 

55. 

-2.5027UC 

60. 

-2.50098C 

65. 

-2.49921E 

70. 

-2.49744E 

75. 

-2.49367E 

80. 

-2.49390C 

85. 

-2.492HE 

90. 

-2.49037C 

95. 

-2.488606 







latitude 


RTMOSPHERIC DENSITY 


MOLECULftR NITROGEN J85.0 KM 10.7 CM FLUX= 92.0’ 


DRY 80. 


I I I i I 

- 

/ H- 80 — -ET \ 



I I I I I ! I I 

,10 .60 .60 t.OO 1.20 1.10 1.60 1.80 

LOCAL TIME 


r Qk 


LEOELS 


SYNBOL 


LN density 


0. 

-2.B2l79C.0i 

5. 

-2.81825C»0I 

10. 

-2.8l470C*0l 

15. 

-2.81116E.01 

20. 

-2.80761E*01 

25. 

-2.80407E.01 

30. 

-2.80052E.01 

35. 

-2.79698C.0l 

40. 

- 2.7934 3C. 01 

45. 

-2.78988C.0l 

50. 

-2.78634E.01 

55. 

-2.782796.01 

60. 

-2.77925C.0l 

65. 

-2.77570C»0l 

70. 

-2.772166*01 

75. 

-2.76861C.0l 

80. 

-2.76507E *01 

85. 

-2.76I52C.01 

90. 

-2.75797E.0l 

95. 

-2.75443C.01 

100. 

-2.75088C.0t 




latitude 


ATMOSPHERIC DENSITY 

MOLECULAR NITROGEN 260.0 Kfl 10.7 CM FLUX= 92.0 DAY 80. 


CONTOUR LEUElS 


to 

.60 


.30 


0.00 


.30 


-.80 



I I I I I I I I I I I 

.20 .40 .60 .80 1.00 1.20 1.40 1.60 1.80 2.00 2.20 10 


SYNeOL LN density 


0. 

- 3.2008 3t *01 

5. 

-3.19370E*01 

10. 

-3.18657E*01 

15. 

-3.17944C .01 

20. 

-3.t723l£»0l 

25. 

-3.165l8E*0l 

30. 

-3.15805C»01 

35. 

-3.15092E.01 

40. 

-3.|4379€»0l 

45. 

-3.13666E.0l 

50. 

- 3.1295 3£. 01 

55. 

-3.12240E.0l 

60. 

•3.11527C.01 

65. 

-3.10814E.01 

70. 

-3.10101E.01 

75. 

-3.09388E.0l 

80. 

-3.086756*01 

85. 

■3.07962E *01 

90; 

- 3.072-4 9£. 01 

95. 

-3.06536C.01 

100. 

-3.058236*01 


LOCAL TIME 




latitude 


PTMOSPHERIC DENSITY 


MOLECULAR NITROGEN 365.0 KM 10.7 CM ELUX= 92.0 


DAY 80. 


CONTOUR LEUELS 


I I 









i I I I I I I I 

.40 .60 .80 1.00 1.20 1.40 1.60 1.80 

LOCAL TIME 


SYMBOL 



LN density 
- 3.640€8E»01 
-3,62850e»01 
-3.61633C*0t 
-3.60415E»0I 
-3.59l97t»0l 
-3.57979€*01 
-3.56761E*0l 
- 3.5554 3E *01 
-3.54325E*01 
-3.53l07E»0l 
-3.51889E*01 
-3.50671E*01 
-3.49453£*0l 
-3.48235e*0l 
-3.47017E*01 
-3.45799E*0l 
-3.4^58lC*0l 
3,433631 *01 
•3.42146€*0l 
-3.40928€*01 
3.397l0£*0l 







ITUDE 






ATnOSPHERIC DENSITY 


CONTOUR UutlS 


MOLECULftR oxygen 


143.0 Kfl 


10.7 CM FLUX= 92.0 


OAT 80. 



SYMBOL lN OtNSITY 


0. 

2.693511 *01 

5. 

-2.691061 *01 

10. 

2.600611 *01 

15. 

-2.686151 *01 

20. 

-2.683701 »01 

, 25. 

-2.601251 *01 

30. 

- 2.670801 *01 

35. 

-2.676341 .01 

40. 

2.673091 .01 

45. 

-2.671441 .01 

50. 

2.660991 .01 

55. 

-2.666531 .01 

60. 

-2.664001 .01 

65. 

2.661631 »0l 

70. 

-2.659181.01 

75. 

-2.656721 .01 

o 

CD 

■2.654271 .01 

05.. 

?.65102( .01 

90." 

2.649371 .01 

95. 

■2.646911 .01 

100: 

2.644461 »01 


LOCAL TIME 



latitude 







LATITUDE 


LOCAL TIME 


0 . 

- 3.459 

5 . 

- 3.450 

10 . 

- 3 . 442 ! 

15 . 

- 3 . 433 ; 

20 . 

- 3.424 

25 . 

- 3 . 416 : 

30 . 

- 3 . 402 : 

35 . 

“ 3.399 

40 . 

- 3 . 390 ! 

45 . 

- 3 . 361 ! 

50 . 

- 3 . 323 : 

55 . 

- 3.364 

60 . 

- 3.3561 

65 . 

- 3 . 342 ; 

20 . 

- 3 . 338 ! 

25 . 

- 3 . 330 ' 

60 . 

- 3.3212 

85 . 

- 3.3131 

90 . 

- 3.3045 

95 . 

- 3 . 295 ! 






A 


rneoL 

LN 0€H 

0. 

-3.96356E 

5. 

- 3.949156 

to. 

-3.93423C1 

15. 

-3.920326* 

20. 

-3.90590e* 

25. 

-3.891496' 

30. 

-3.822086' 

35. 

-3.862666 

40. 

-3.048256 

45. 

-3.833836 

50. 

- 3.8194 2E* 

55. 

- 3.805016 ♦ 

60. 

-3.2^9E( 

65. 

-3.226186 < 

20. 

-3.261266* 

25. 

-3.242356 

80. . 

-3.232946' 

85. 

-3.218526* 

90. 

- 3.204116 •( 

95. 

-3.689696' 





LATITUDE 




ftTOmC OXYGEN 

I I 

5 


ATMOSPHERIC DENSITY 

119.0 KM 10.7 Cn FLUX* 92,0 

till 

5 S ~ 5 - 



I I I I t I 

1.20 1.40 1. 

LOCPIL TIME 


,40 


.60 


.80 


1.00 


DAY 00. 


tilt 



CONTOUR LEuCLS 


SYMBOL LN OCNSITt 


0 . 

'2.72409E*01 

5. 

-2,72291E*01 

10 . 

-2.72173E»01 

15. 

'2.72055£*0l 

20 , 

-2.7l937t*0l 

25. 

-2.7l0l8E*Ot 

30. 

'2.71700e*0l 

35. 

-2.71502E*OI 

40. 

-2.71464E»01 

45. 

-2.71346E»01 

50. 

-2.71228E.01 

55. 

-2.71109C*01 

60. 

-2.7099tE*0l 

65. 

-2.70873E.01 

70. 

-2.70755C.0l 

75. 

-2.70G37E.Ot 

BO. 

-2.70518C.01 

85. 

-2.70400E»01 

90, 

-2.70282C.01 

95. 

-2.70164C .61 

100 . 

-2.70046C.01 




latitude 






latitude 





riTUDE 


ERIC DE 





LATITUDE 


ATMOSPHERIC DENSITY 


CONTCXJR LEUCLS 



ATOMIC OXYGEN 365.0 Ktl tO,7 CM FLUX= 92.0 DAY 80. 


I I I I t I I I I I 



LOCAL TIME 


SYMBOL LN OEMS IT Y 
0. -3,5t;^92E*01 

5. -3.5J28^*OJ 

10. -3.50;y8i£,oI 

15. -3.50276£*01 

20. -3.49;Y7ie*oi 

25. -3.49265€»01 

30. -3.48260C*0I 

35. -3.40255t*Ol 

40. -3.42;»5OC*0l 

45. -3.42244E*0l 

50. “3.46239C*01 

55. -3.46234E»0l 

60. -3.45229€»0l 

65. -3.45223C»01 

PO. -3.447I8£*0J 
25. -3.442l3E.0l 

90. -3.43702E*0l 

95. -3.43202C.01 

90. -3.42692E.0l 

95. -3.42192E.0l 

100. -3.4te86E.01 


'£) hrl 
<T> 

Oq 



latitude 


PHERIC DENSITY 

KM 10.7 CM FLUX= 92.0 


DOT 60. 

I 


CONTOUR LEUCLS 


0. 

-3.78118€< 

5. 

-3.77349C 

10. 

- 3,765806 

15. 

-3.75811E* 

20. 

-3.750416 

25. 

-3.74272E 

30. 

-3.735036 

35. 

-3.727346 

^0. 

-3.719656 

45. 

-3.711966' 

50. 

-3.704276 

55. 

-3.696576 

60. 

-3.688886 

65. 

- 3.681196 < 

70. 

-3.673506 

75. 

-3.665816 

80. 

- 3.658126 < 

85. 

-3.65042E 

90. 

-3.642736 

95. 

-3.63504E 






LATITUDE 







RinOSPHERIC DENSITY 

nOLECULAR NITROGEN H3.0 KN 10.? CM FLUX= 92.0 


DAY 172. 


CONTOUR levels 



1 I I I I ) 

.40 XO .80 iUM 1.20 1.40 

LOCAL TIME 


SYflGOL 
. 0 . 

5. 

10 . 

15. 

20 . 

25. 

30. 

35. 


2M 


LN OENSIIT 
-2.51902E*01 
-2.51743E.01 
-2.5158 3E»0l 
-2.5l424E*0l 
-2.51264E*01 
-2.5ll04E*01 
-2.50945£*0l 
-2.50785E»0I 
-2.50626£*01 
-2.50466E*01 
-2.50307E»0l 
-2.50147E.01 
-2.49988E.01 
-2.49828E.0l 
-2.49668E.01 
-2.49509E.0l 
-2.49349E*01 
-2.49190E .01 
-2.49030C.01 
-2.48871E.01 
-2.48711E.01 





ATMOSPHERIC 


flOLECULflR NITROGEN 185.0 Kfl 10.7 CM I 


i I I I I 






latitude 


CONTOUR L,Eu€lS 







ninOSPHERIC DENSITY 

HOLECULAR NITROGEN 365.0 KM 10.7 CM FLUX= 92.0 

I I I I I t 




DAY 172. 



I i i I I I 

,40 .60 .00 1.00 - 1.20 1.40 

LOCRL TIME 


CONTOUR LEUCLS 

SYMBOL LN density 
0. -3.63692E*01 

5. -3,62534E*01 

10. -3.6137;'E*01 

15. -3,60219E«01 

20. '3.59062E»01 

25. -3.57905E*01 

30. -3,567V€»01 

35. -3.55590E»01 

40. ■ '3.54433E»01 
45. -3.53275C*01 

50. -3.52118E*01 

55. -3.50961E*0l 

60. -3.43803E*01 

65. -3.486466*01 

70. -3.474096*01 

75. -3.46331E*0l 

80. -3.45174E*01 

85. -3.440176. 0) 

90. 3.42859C *01 

95. -3.417021*01 

100. -3.405446 * 01 




latitude 




LATITUDE 






latitude 


ATMOSPHERIC DENSITY 


, don T OUR' tEUCLS 


MOLECULAR OXYGEN 

I I 


185.0 KM 


10.7 CM FLUX= 


DAT 172. 

I 






I I I 

1.00 1.20 U40 

LOCAL TlhE 


SYMBOL 


LN OCNSI-T t 



0. 

- 3;01601E'»01 

5. 

- 3.012 36E *01 

10. 

-3.OO07OE*Ol 

15. 

-3.00505E»01 

20. 

-3.00139E*01 

25. 

-2.99774E»01 

30. 

-2.99408E*01 

35. 

-2.99043E»0l 

40. 

-2.98677£*0l 

45. 

-2.9§'3t2E*0l 

50. 

-2.97946E»0l 

55. 

-2.97581E»01 

60. 

-2.972156 »0l 

65. 

-2.96850£»0l 

70. 

-2.96404E»Ol 

75. 

-2.96ll9e»0l 

80. 

-2.957536*01 

85. 

•2,95388t*0l 

90. 

-2.950236*01 

95. 

-2,946576*01 

too. 

-2.94292E»01 





latitude 


ftinOSPHERIC DENSITY 

molecular oxygen 260.0 KM 10.7 CM FLUX= 32.0 


DAY 172. 


CONTOUR lEvClS 



.40 .60 .80 1.00 1.20 


SYMBOL LN DENSlT’f 
0. -3.4460lE*Ol 

5. -3.43324f0l 

10. -3.43166E»0l 

15. -3.42409E*01 

20. -3.4l652E*0l 

25. -3.40894E*01 

30. '3.40l37E*0t 

35. -3.39373E*0l 

40. -3.38622E»01 

45. -3.37865E»01 

50. -3.37107E»01 

55. -3.36350E«01 

60. -3.35593E»0l 

65. -3.34835E»01 

70. -3.34078E*01 

75. -3.33320£»01 

80. -3.32563E»oT" 

85. -3.31806E*01 

90. -3.31048C»01 

95. -3.30291E»01 

100. -3.29534E»01 


2.20 to 


LOCAL TIME 




TUDE 





latitude 


CONTCXJR LE«JClS 















latitude 





latitude 





ninOSPHERIC DENSITY 


C0NT0O*f lEWClS 



.40 .60 .60 1.00 1.20 1.40 1.60 1.80 2.00 2.20 10 

LOCAL TIflE 


SYMBOL LM DCNSITt 


0. 

•3.28l23E‘0t 

5. 

-3.2.»252t*0l 

10. 

•3.2;>38lt*01 

15.' 

-3.27009t*01 

20. 

-3.26638E«01 

25. 

-3.2626^E*0l 

30. 

-3.25895E»01 

35. 

-3.25524E*0l 

40. 

-3.25152E»01 

45. 

-3.2478IE»01 

50. 

'3.24410€*01 

55. 

-3.?403BE»0l 

60. 

- 3.2366;>E *01 

65. 

- 3.2329€f *01 

^0. 

3.22924E *01 


3.22553E *01 

80. 

3.22182E *01 

85. 

3.21810E*01 

90. 

-3.21439C.oi 

95. 

■ 3.210S;’C *01 

100. 

■ 3.20696E *01 








latitude 



E 



LM D€^ 

0. 

-3.29407E 

5, 

-3.;'0491E 

10. 


15. 

-3.;'G660e< 

20. 

-3.25744E 

25. 

- 3.74829C 

30. 

- 3.7391 3£< 

35. 

-3.72997E 

<0. 

■3.72082C 

45. 

-3.71J66E* 

50. 

- 3.70250C • 

55. 

- 3.69335E 

60. 

- 3.604 19E. 

65. 

-3.67503E< 

;^o. 

- 3.66588E 

/5. 

-3.65672C 

80. 

- 3.64756r 

85. 

-3.63041E* 

90. 

- 3.62925C • 

95. 

• 3.62009E • 

JO. 

3.61094E *0 





latitude 




LATITUDE 







LATITUDE 


SPHERIC DENSITY 

.0 KM 10.7 CM FLUX= 92.0 

I t 


DAY 172. 


CONTOUR LEUeCS 


.30 — 


0. 

-3.20197E^ 

5. 

-3.19430E 

10. ' 

-3.18679e< 

15. 

-3.17919E + 

20. 

-3.17160£» 

25. 

-3.16400E^ 

30. 

-3.1564 IE < 

35. 

-3.14882E 

^0. 

-3.14122E* 

45. 

-3.1336 3£ 

50. 

-3.12604E4 

55. 

-3.11844E< 

60. 

-3.1108SE* 

65. 

-3.10325E^ 

70. 

-3.09566E 

75. 

-3.08807E' 

eo. 

-3.08047C< 

85. ■ 

"3.07288E 

90. 

-3.06529E- 

95. 

-3.05769C 






PiTnOSPHERIC DENSITY 


MOLECULAR NITROGEN 


365.0 KN 


10.7 CM FLUX= 92.0 


DAY 172, 







lptitude 


CONTOUR LEUCLS 










LATITUDE 


ATMOSPHERIC DENSITY 


MO.ECULAR OXYGEN 


119.0 KM 10.7 CM FLUX= 92,0 


DAY 172. 


CONTOUR LEUELS 



SYMBOL LN DENSITY 
0. -2.^S266E-»01 

5. -2.45085E*0I 

10. -2.44903E*01 

15. -2.44722E*0l 

20. -2.44540E>01 

25. -2.44359E*01 

30. -Z.4^17BE»01 

35. -2.43996E+01 

40. -2.43815E+01 

45. -2.43833E+01 

50. -2.43452E*0l 

55. -2.43271E*01 

60. -2.43089£»0l 

65. -2.429O0E*O1 

70. -2.42726E+01 

75. -2.42545E*0l 

80. -2.42364E+01 

85. -2.42182E*01 

90. -2.42001E*01 

95. -2.41820£*01 

100. -2.4l638C»0l 


.40 .60 


.80 1.00 1.20 1.40 1.60 1.80 

LOCAL TIME 





latitude 






ATnOSPHERIC DENSITY 


MOLECULAR OXYGEN 


185.0 KM 


10.7 CM FLUX= 92.0 


□AY 172, 





LATITUDE 












LATITUDE 


OTOfllC OXTCEN 


ATMOSPHERIC DENSITY 


119.0 KM 10.7 CM FLUX= 92.0 


DAT 172. 


5 W 



CONTOUR LCUEUS 


SYMBOL LN density 
■ 0. -2.73503E*01 

5. -2.73329E*01 

10. -2.73155E+01 

15. -2.72981E+01 

20. -2.72807E»01 

25. -2.72632E?0l 

30. -2.72458E*0l 

35. -2.7228-^E*0l 

■♦0. -2.72110E+01 

^5, -2.71936E+01 

50. -2.71762E*01 

55. -2.7l58BE»0l 

60. -2.7H13E*01 

65. -2.71233E*01 

70. -2.710S5E*01 

75. -2.7069lE*01 

80. -2.70717E*01 

85. -2.70543E»01 

90. -2.70368E+01 

95. - 2.701 94E *01 

100. -2.70020E«C1 


.40 .60 


LOCAL TIME 




latitude 


ATOMIC OXYGEN 


RTdOSPHERIC DENSITY 


H3.0 KM 10.7 CM FLUX« 92.0 


DAY 172. 


CONTOUR LEUELS 



I I 

.40 .60 


SYMBOL LN density 
0. -2.89543E*0l 

5. -2.89286E*01 

10. -2.89028E«01 

15. -2.88770E*01 

20. -2.88512E*01 

25. -2.88255£*01 

30. -2.87997E+01 

35. -2.87739E*01 

40. -2.87481E*01 

45. -2.87224E*0l 

50. -2.86966E«0I 

55. -2.86708E*01 

60. -2.86451E*01 

65. -2.86193E*01 

70. -2.85935E«01 

75. .-2.85677E»0l 

80. -2.85420E*01 

85. -2.85162E*01 

90. ■2.84904E«01 

95. -2.84647E*0l 

100. -2.84389E*01 


2.20 10 


LOCAL TIME 





RTMOSPHERIC DENSITY 


ATOMIC OXYGEN 


185.0 KM 


10.7 CM FLUX* 92.0 


DAY 172, 






LATITUDE 







latitude 







latitude 







LATITUDE 









latitude, 

^ « 






ATMOSPHERIC DENSITY 


nOLECULAR NITROGEN 


3S5.0 KN 


10.7 CM FLUX= MO.O 


OAT 172, 




latitude 


ERIC DE 



rnSOL 

LN D£NSI 

0. 

-3.83222E*0l 

5. 

-3.81803E*01 

10. 

“3.8O303£*O1 

J5. 

-3.2896-1E*01 

20. 

-3.225-14E*0l 

25. 

-3.26l25E*0l 

30. 

- 3.2-1 206E* 01 

35. 

-3.23286E*01 

■<0. 

-3.21862E*0J 

-15. 

-3.20i42E»01 

SO. 

-3.69028E*0J 

55. 

-3.62609E*01 

60. 

-3.66109E*O1 

65. 

-3.6-1220E*Ol 

20. 

-3.63350£*0J 

25. 

-3.6193lE*0l 

00. 

-3.605l_2E»01 

85. 

-3.59092E»01 

90. 

- 3.5262 3£ *01 

95. 

-3.56253£*0l 





ATMOSPHERIC DENSITY 


MOLECULAR OXYGEN 


119.0 <M 


10.7 CM FLUX= 


day 172. 


CONTOUR LEUELS 



SYMBOL LN density 
0, 'S,979OOE*02 

5. -5.97882E*02 

10. -5.978G4E»02 

15. -5.978^GE*02 

20. -5.97828E»02 

25. - 5.978 lOE *02 

30. -5.97792E*02 

35. -5.97?73£*02 

-TO. -5.97755E»02 
45. -5.97737E*02 

50. “5.97718E»02 

55. -5.97701E*02 

60. -5.97e83E*02 

65. -5.97665E*02 

70. -5.97646E*02 

75. '5.9762BE*02 

80. - 5.976 lOE *02 

85. -5.97592E »02 

90. -5.9757-»E»02 

95. -5.97556E»02 

100. - 5.975 37E* 02 


.40 .60 


LOCAL TIME 






latitude 


flinOSPHERIC DENSITY 


MOLECULAR OXYGEN 


143.0 KM 10.7 CM FLUX= 140.0 


DAY 172. 


CONTOUR LEuELS 



SYT«OL LN density 
0. -7.40472E*02 

5. -7.40442E*02 

10. -7.40412E*02 

15. -7.40382E»02 

20. -7.40352E»02 

25. -7.40322E*02 

30. -7.40292E*02 

35. -7.40263E*02 

40. -7.40233E*02 

45. -7.40203E *02 

50. -7.40173E*02 

55. -7.40143E*02 

60. -7.40113E»02 

65. -7.40083E*02 

70. -7.4005 3E*02 

75. -7.40023E*02 

80. ■ -7.39993£»02 
85. -7.39363E*02 

90. -7.39933E*02 

95. - 7.3990 3£ *02 

100. -7.39873E*02 


.40 .60 .80 1.00 1.20 1.40 

LOCAL TIME 







latitude 









latitude 






ftTOmC OXTCEN 


RTflOSPHERIC DENSITY 


119.0 KM 10.7 CM FLUX= HO.O 


DOT 172. 


CONTOUR LEUELS 



StMBOL 


75 _ 


LN DENSITT 
-2.73359E*0l 
-2.73166E»01 
-2.72973E*01 
-2.72781E*0l 
-2.72580E*O1 
-2.7239SE»01 
-2.72203E*0l 
-2.72011E*OI 
-2.71818E*01 
-2.7162eE»01 
-2.7M33E«01 
-2.712^1E'*»0l 
-2.71048E*01 
-2.70856E»01 
- 2.7066 3£»01 
-2.7047lE»qi 
-2.70278E«01 ' 
-2.7OO06E»Ol 
-2.69893E *01 
-2.69700E»01 
-2.69508E*01 


LOCfiL TIME 





latitude 


ftTOHIC OXYGEN 


AinOSPHERIC DENSITY 

H3.0 KM 10.7 CM FLUX= HO.O 

I I I I 


DRY 172. 



CONTOUR LEUCLS 

SYMBOL LN density 
0. -2.80992E»Ol 

5. -2.80713E*O1 

10. -2.8843-IE^Ol 

15. -2.8815-1E*0l 

20. -2.87075E*Ol 

25. -2.87596E*01 

30. -2.87316E*01 

35. -2.87037E*01 

“^O. -2.06758E»O1 

45. -2.06470E*O1 

50. -2.06199E*O1 

55. -2.05919E*Ol 

60. -2.0564O€*O1 

65. '2.05361E*OI 

70. -2.85081E*01 

,75. -2.040O2E*O1 

00. -2.04523E»Ol 

85. -2.04243E*O1 

90. -2.03964E»O1 

95. -2.03805E*O1 

100. -2.83405E»01 


1.00 1.20 1.40 

LOCAL TIME 


1.60 1.80 







ATOMIC oxygen 


RTMOSPHERIC DENSITY 


185.0 KM 10.7 Cfl FLUX= HO.O 


DRY 172. 


CONTOUR LEUELS 



SYMBOi 

0 . 

5. 


LN DENSIT 
•3.07237E.01 
-3.0G907E*01 
-3,06578E»01 
-3.062^8E»0l 
-3.05918E*01 
-3.05588£*0l 
-3.05259E*01 
-3.0i929C»0t 
-3.0^599E‘01 
-3.0'<269E*01 
-3.03940E*01 
-3.036l0C*0l 
-3.O320OE*Ol 
-3.02950E*01 
-3.02620e*01 
-3.02291E*01 
-3.01961E*01 
-3.0lB31E»0l 
-3.0t301E*01 
-3.00972E*01 
-3.00642E«01 


.40 .60 


t .20 1.40 

LOCAL TIME 


1.60 1.80 




latitude 






pthospheric density 


fiTomc OXTCEN 


365.0 KM 


10.7 CM FLUX* HO.O 


DAY 172. 









